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■ Abstract 

A search for — > oscihations has been conducted with the LSND ap- 
^ \ paratus at the Los Alamos Meson Physics Facihty. Using from vr^ de- 

5h \ cay in flight, the v^, appearance is detected via the charged-current reaction 

f e C — > e~ X. Two independent analyses observe a total of 40 beam-on 
high-energy electron events (60 < Eg < 200 MeV) consistent with the above 
signature. This number is significantly above the 21.9 it 2.1 events expected 
from the fe contamination in the beam and the beam-off background. If 
interpreted as an oscillation signal, the observed oscillation probability of 
(2.6 it 1.0 lb 0.5) X 10"'^ is consistent with the previously reported Ue 
oscillation evidence from LSND. 

14.60.Pq, 13.15.-bg 
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In this letter we describe the results of a search for —>■ p^. oscillations using a flux 
from 7r+ decay in flight (DIF). The data were taken with the Liquid Scintillator Neutrino 
Detector (LSND) at the Los Alamos Meson Physics Facility (LAMPF). The result of a 
search for v^. oscillations, using a flux from /i"^ decay at rest (DAR), has already 

been reported in Ref. where an excess of events was interpreted as evidence for neutrino 
oscillations. The analysis presented here uses a different component of the neutrino beam, a 
different detection process, and has different backgrounds and systematics from the previous 
DAR result, providing an independent check on the existence of neutrino oscillations. 

The primary source of DIF for this experiment is the A6 water target of the LAMPF 
800 MeV proton linear accelerator. Approximately 3.4% of the 7r+ produced in the 30-cm 
long target decay in flight before reaching the water-cooled copper beam stop, situated 1.5 m 
downstream. The generated flux, with energies up to 300 MeV, is illustrated in Fig. |l|(a), 
as calculated at the center of the detector, 30 m away from the beam stop. Two upstream 
thin carbon targets, Al and A2, located at 135 m and 110 m from the detector center, 
respectively, provide additional small contributions to the Vy, flux - also shown in Fig. |l](a). 
However, for z/^ v^, oscillations with small Am^ values, the z/^ flux from Al and A2 can 
have a significant effect due to the longer baselines. The main beam-related backgrounds 
(BRB) to the search come from the intrinsic u^. component of the beam, shown in 

Figs. |l](b) and (c). The flux from tt^ e^Ve DIF is suppressed by the branching ratio of 
1.24 X 10~^, while the flux from /i+ e^z^gZ/^ DIF is suppressed by the longer /i lifetime 
and the kinematics of the three-body decay. The neutrino flux calculations are described 
in detail in Ref. and yield a systematic error of 15% for the DIF flux. The LSND 
measurement of the exclusive reaction z/^ — > /x^ ^"^Ngs, which has a well understood cross 
section, confirms the calculated flux to within a 15% statistical error 0. 

The data discussed here correspond to 14772 Coulombs of protons on target (POT) 
during the years 1993 (1787 C), 1994 (5904 C), and 1995 (7081 C). The beam duty factor - 
defined as the ratio of data collected with beam on to that with beam off - has a weighted 
average of 0.07 for the three years. 

The LSND apparatus, described in detail elsewhere f^, consists of a steel tank filled with 
167 metric tons of liquid scintillator and viewed by 1220 uniformly spaced 8" Hamamatsu 
photomultiplier tubes (PMT). The scintillator medium consists of mineral oil {CH2) with 
a small admixture (0.031 g/1) of butyl-PBD. This mixture allows the detection of both 
Cerenkov and isotropic scintillation light, so that the on-line reconstruction software provides 
robust particle identification (PID) for electrons, along with the event vertex and direction. 
The electronics and data acquisition (DAQ) systems were designed to detect and record 
related events separated in time. 

Despite 2.0 kg/cm^ shielding above the detector tunnel, there remains a very large back- 
ground to the oscillation search due to cosmic rays. This background is highly suppressed 
by a veto shield 0, which provides both passive and active shielding. It is viewed by 292 
uniformly spaced 5" EMI PMTs and has a threshold of 6 PMT hits. Above this value a 
signal holds off the trigger for 15.2 /is while inducing an 18% dead-time in the DAQ. A veto 
inefficiency < 10~^ is achieved off-line with this detector for incident charged particles. The 
veto inefficiency is obviously much larger for incident cosmic-ray neutrons. 

A GEANT-based Monte Carlo (MC) is employed to simulate interactions in the LSND 
tank and the response of the detector system. The detector response parameterizations 
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were measured either in a test beam or in a controlled setting. The electron simulation 
is calibrated below 52.8 MeV using Michel electrons from the decay of stopped cosmic-ray 
muons and then extrapolated into the DIF energy range. The MC data set used to calculate 
electron selection efficiencies in the DIF analysis (DIF-MC) uses the calculated flux, 100% 
~^ transmutation, and the Vf,C e' X cross section calculated in the CRPA model 

i- 

Candidate events for oscillation from the DIF z/^ flux consist of a single, isolated 

electron (from the Vf.C X reaction) in the energy range 60-200 MeV. The lower limit 
is chosen to be well above the endpoint of the Michel electron spectrum (52.8 MeV) to avoid 
backgrounds induced by cosmic-ray muons and beam-related and events. The upper 
limit of 200 MeV is the energy above which the beam-off background rates increase, and the 
expected signal becomes much attenuated. 

A preliminary selection was made to arrive at an initial data sample. The electron 
PID parameters used in the DAR analysis retain high efficiency (98.1±1.7%), but have 
limited background rejection in the DIF energy range. New PID parameters developed for 
this analysis are used in final event selection as described below. To reduce the cosmic-ray 
muon related background several cuts are made. First the veto shield is required to have 
less than 4 active PMTs. Second, the events must be reconstructed more than 35 cm from 
the surface defined by the inner PMT faces. Third, events with another event of 200-700 
hit PMTs in the following 30 //s (characteristic for Michel electrons) are dropped if they 
are within 200 cm of the subsequent event, or have more than 600 hit PMTs, eliminating 
stopping muons not vetoed by the DAQ. Finally, the event must not have had any previous 
activity above 600 hit PMTs or within 200 cm in space in the preceeding 30 /is, removing 
remaining Michel electrons. These cuts have an overall efficiency of 82.4±2.7% for electrons 
in the DIF energy range. 

The event reconstruction and PID techniques used in the DIF analysis were developed 
to utilize fully the capabilities of the LSND apparatus. The basis for the reconstruction 
is a simple single track event model, parametrized by the track starting position and time, 
direction, and energy. For any given event, the expected photon intensities and arrival times 
are calculated from these parameters at all PMTs. A likelihood function that relates the 
measured PMT charge and time values to the predicted values is used to determine the best 
possible event parameters and also provides PID. Two independent analyses were performed, 
sharing basic goals, but differing in approach and parameterizations. Both analyses are 
described in detail elsewhere 0. 

The essential goal of both analyses is to select events consistent with DIF candidate 
electrons, while eliminating remaining backgrounds from cosmic-ray interactions, including 
neutrons and photons. The electron identification relies primarily on the differences in the 
timing characteristics of the components of light produced in the event: scintillation light, 
and Cerenkov light, both direct and rescattered - Figs, ^(a) and (b). Furthermore, the 
event likelihood fitting returns also the fraction of direct Cerenkov light in the event, which 
provides excellent rejection against neutrons - Fig. |^(c). 

High energy 7 rays, from 7r° produced by neutron interactions in the lead shielding of 
the veto system, enter the detector fiducial volume without leaving a veto signal. The 
charged particles resulting from their interactions in the liquid point predominantly into 
the detector volume and are difficult to distinguish from electrons from the Vf.C —>■ e~ X 
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reaction on the basis of electron PID alone. The backwards projected track-length to the 
edge of the detector volume, s, is used to remove these events, concentrated at low values 
of s - Fig. ^(d). Events with any veto hits in time with the event, and along the backward 
extrapolation of the track, are also rejected. 

Finally, the electron events in the final DIF sample are required to have cos Of. < 0.8, 
where de is the angle between the reconstructed direction and the incident neutrino beam. 
This greatly reduces the BRB from the forward-peaked t'^e elastic scattering, while retaining 
a high efficiency for the DIF signal. 

After applying all of the respective selection criteria, both analyses obtain a significant 
and consistent beam-related event excess. One analysis ends up with 23 beam-on events 
and 114 beam-off events (8.0 rescaled for the duty factor), which corresponds to 15.0 excess 
events. The other analysis ends up with 25 beam-on events and 92 beam-off events (6.4 
rescaled), which corresponds to 18.6 excess events. Their efficiencies are 8.4% and 13.8%, 
respectively, calculated for to the d > fiducial volume. 

As already mentioned, the main BRBs in the DIF oscillation search come from the 
intrinsic Uf. contamination in the beam, /i"*" e'^Uf.iy^ (DIF), and e^z/g (DIF). These 

backgrounds are calculated using the beam MC neutrino fiuxes and the u^C cross section 
calculated in the CRPA model. The z/^e elastic scattering background from the DIF fiux 
is greatly reduced by requiring cos^e < 0.8. The last relevant background, 7r+ IJ-^^fi 
DIF followed by ly^C i^f^Cn" coherent scattering, is calculated using the cross section in 
Ref. 0]. Backgrounds from the f^C fi^X reaction are negligible. The four relevant BRBs 
are summarized in Table |. The total BRBs calculated for the two analyses yield 4.5 and 8.5 
events, respectively, which thus leaves a significant excess of events (10.5 and 10.1 events, 
respectively) above the expectation from conventional processes. The probabilities that the 
number of expected background events (12.5/14.9) fiuctuate up to the observed beam-on 
numbers (23/25) are 0.7 x 10^^ and 1.6 x 10~^, respectively. 

Since both analyses have low efficiencies, different reconstruction software, and different 
selection criteria, the two samples need not necessarily be identical. Both the logical AND 
and OR of the two samples have been extensively studied in MC simulations and the results 
are consistent with the expectations. For the final DIF sample we have elected to use the 
logical OR of the events. This procedure minimizes the sensitivity of the measurement to 
uncertainties in the efficiency calculations, is less sensitive to statistical fiuctuations, and 
also yields a larger overall efficiency. Table |I| summarizes the final event samples for the 
individual analyses, for their overlap (AND), and for the final sample (OR). 

In the following we interpret the observed event excess of the OR sample in terms of the 
simplest, two-generation mixing neutrino oscillations formalism. In this model the oscillation 
probability is given by 

P = sin^ 26 sin^ (^1.27 ^m^^^ , (1) 

where 9 is the mixing angle, Am^ (^V^/ c^) is the difference of the squares of the masses of the 
appropriate mass eigenstates, L (m) is the distance from neutrino production to detection, 
and El, (MeV) is the neutrino energy. Since the distance to the source is ambiguous because 
of the presence of multiple beam targets (Al, A2, and A6), the energy distribution alone is 
used to determine the confidence levels in the (sin^ 29, Am^) parameter space. Fig. ^ shows 
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the 95% confidence level contours that result from the fit. This result is consistent with the 
previous LSND DAR result shown superimposed in Fig. 0. The oscillation probability 
is (2.6 ± 1.0 ± 0.5) X 10~^, where the second error is systematic, as described below. 

The neutrino cross sections and fluxes constitute the largest source of systematic uncer- 
tainty for the DIF analysis. Although our measurement of V(.C scattering using the DAR 
z/g flux agrees well with calculations 0, our measurement of the inclusive v^C cross section 
is 45% below the CRPA calculation The flux in the f^C measurement is the same 
as for this DIF oscillation analysis, and it is possible that the V(.C cross section at these 
higher energies also is below the CRPA calculation. The expected average number of events 
is given by 

Ntotai = e a^^c {^v^Pv^^u, + ^u,) + Nbub (2) 

where e is the event selection efficiency, ^ue/u^ are the neutrino fluxes, a^^c is the neutrino 
cross section, and Nbub is the beam-unrelated background (BUB). The oscillation signal 
is proportional to the same product {e a^^c ^u^) as the neutrino background, since is 
proportional to . The effect of lowering the product e $ a^^c is to reduce the predicted 
BRB, which raises the observed oscillation signal. Only by raising the product e $ a^^c is 
the oscillation signal decreased. In order to calculate conservative confidence regions, we 
allow the value of e^a^^c to vary between 21% above to 45% below the calculated value. 
Only a symetrical 21% systematic error is used in the oscillation probability. 

We have described a search for z/g C e~ X interactions for electron energies 60 < 
Ef, < 200 MeV. Two independent analyses observe a number of beam-on events significantly 
above the expected number from the sum of conventional beam-related processes and cosmic- 
ray (beam-off) events. The probability that the 21.9 ± 2.1 estimated background events 
fluctuate into 40 observed events is 1.1 x 10"'^. The excess events are consistent with z/^ — >■ 
Ue oscillations with an oscillation probability of (2.6 ± 1.0 ± 0.5) x 10~^. A fit to the 
energy distribution events, assuming neutrino oscillations as the source of /^e, yields the 
allowed region in the (sin^ 29, /S.m?) parameter space shown in Fig. |^. This allowed region 
is consistent with the allowed region from the DAR search reported earlier by LSND. This 
DIF oscillation search has completely different backgrounds and systematic errors 
from the — > Ve DAR oscillation search and provides additional evidence that both effects 
are due to neutrino oscillations. 
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TABLES 

TABLE I. Background estimates for the oscillation search for the d > fiducial 

volume, 9.2 x 10^^ POT, and for electron energies between 60 MeV and 200 MeV. These numbers 
are illustrative for an electron selection efficiency of 0.10, independent of energy. The actual 
efficiencies in the two analyses are slightly different and energy dependent. 



Process 




Flux 


< 0- >^ 


Eff. 


Number 






(cm-2pOT-i) 


(10-40 cm-2) 


(%) 


of Events 


VeC — > eT 


X ill DIF) 


3.8 X 10-^4 


28.3 


10.0 


3.8 


VeC — > e~ 


-X (tt DIF) 


8.3 X lO"^-"^ 


79.2 


10.0 


1.6 


v^C ^ v^Ct^° 


6.5 X 10"^^ 


1.6 


6.0 


0.3 


v^ 


e 


6.5 X 10"^^ 


0.00136 


0.5 


0.1 


Total background 








5.8 


TABLE II. Comparison of results for the two analyses (labeled here 


as A and B), their logical 


AND and OR. All errors 


are statistical. 








Data 


Beam 


BUB BRB 


Excess 


Eff. 


Osc. Prob. 


Set 


On/Off 






(%) 


(xlO-3) 


A 


23/114 


8.0 ±0.7 4.5 ±0.9 


10.5 ±4.9 


8.4 


2.9 ± 1.4 


B 


25/ 92 


6.4 ±0.7 8.5 ±1.7 


10.1 ±5.3 


13.8 


1.7 ±0.9 


AND 


8/ 31 


2.2 ±0.3 3.1 ±0.6 


2.7 ±2.9 


5.5 


1.1 ± 1.2 


OR 


40/175 


12.3 ±0.9 9.6 ±1.9 


18.1 ±6.6 


16.5 


2.6 ± 1.0 
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FIG. 1. Calculated and DIF fluxes at the detector center from the A6 target (solid 
histograms) and from the A1+A2 targets (dashed histograms). 
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FIG. 2. Timing likelihoods for (a) the entire event and (b) the Cerenkov region only, (c) 

Cerenkov-to-scintillation density ratio, p. (d) Projected track-length to the tank wall intersection. 
(a)-(c) correspond to all (beam on+ofF) DIF data after the prc-selection and (d) after all other cuts 
have been applied. All superimposed distributions (dashed) correspond to the DIF-MC simulation, 
normalized to the same areas. 
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FIG. 3. The 95% confidence level region for the DIF i/^ —>■ Ug along with the favored regions 
from the LSND Vg DAR measurement (dotted contours). 
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